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Abstract 

The Next-to-Minimal Supersymmetric extension of the Standard Model (NMSSM) fea¬ 
tures extra new sources for CP violation. In contrast to the MSSM CP violation can already 
occur at tree level in the Higgs sector. We investigate the range of possible allowed CP- 
violating phases by taking into account the constraints arising from the measurements of 
the Electric Dipole Moments (EDMs) and the latest LHC Higgs data. Our analysis shows 
that large CP-violating phases, that are NMSSM-specific, are not in conflict with the EDMs. 
They are dominantly constrained by the Higgs data in this case. We use our results to 
investigate the prospects of measuring CP violation through the combined measurement of 
Higgs rates, on the one hand, and in observables based on CP-violating Higgs couplings to 
tau leptons on the other hand. 
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1 Introduction 


The discovery of a Standard Model (SM)-like Higgs boson with mass around 125 GeV by the 
Large Hadron Collider experiments ATLAS and CMS represents a milestone for particle 
physics. This discovery supports the Higgs mechanism which allows for the generation of particle 
masses without violating the underlying gauge symmetries of the SM. It has long been realized 
that the maximum allowed symmetry compatible with space-time symmetry is supersymmetry 
(SUSY) |3j, which relates bosonic and fermionic degrees of freedom. Due to this and many 
other virtues as well, SUSY has become one of the most popular and most intensely studied 
symmetries beyond the SM (BSM). While so far no direct sign of new physics has been found, it 
remains possible that the SM-like Higgs boson discovered at the LHC is in fact a SUSY Higgs. 

The Minimal Supersymmetric Standard Model (MSSM) [4], requires at least two complex 
Higgs doublets. In the Next-to-Minimal SUSY model (NMSSM) |5|| another complex singlet 
superfield is added to the Higgs sector. The coupling of the singlet field to the MSSM Higgs 
doublets allows for a dynamical solution of the n problem |6[ when the neutral component of the 
singlet field acquires its vacuum expectation value (VEV). The totality of ten degrees of freedom 
of the Higgs doublet and singlet fields leads to seven physical Higgs bosons after electroweak 
symmetry breaking. In the CP-conserving NMSSM these are three CP-even and two CP-odd 
neutral Higgs bosons plus two charged ones, whereas in the CP-violating case all Higgs bosons 
mix and do not carry a definite CP quantum number any more. Besides the direct detection 
of more Higgs states, an extended Higgs sector can manifest itself in modified Higgs couplings 
of the SM-like Higgs boson. These arise from the mixture with other Higgs states or from new 
physics effects induced through radiative corrections to the Higgs couplings and/or in the loop 
mediated Higgs interactions with the photons and gluons. Furthermore, Higgs decays into other 
lighter non-SM particles can be realized leading to modified branching ratios, including also the 
possibility of invisible decays. 

CP violation is one of the three Sakharov conditions [T] for baryogenesis, leading to matter- 
antimatter asymmetry in the universe. In the SM the only source of CP violation is given by 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix |8j 91. While the SM provides the necessary 
ingredients for all Sakharov conditions, CP violation based on the CKM matrix is too small to 
explain quantitatively the observed asymmetry^ This motivates studying BSM theories which 
include additional sources of CP violation. SUSY contains many new sources of CP violation. 
In the MSSM CP violation in the Higgs sector itself cannot occur at tree level and is radiatively 
induced. In the NMSSM Higgs sector, however, CP violation can already show up at tree level 
through NMSSM specific complex couplings which induce CP-violating doublet-singlet mixing. 

In this paper we shall be concerned with the CP-violating NMSSM, including CP-violating 
effects in the phenomenology of the 125 GeV Higgs boson as well as other effects in the Higgs 
sector as follows. For example, in the CP-violating NMSSM, a non-zero CP-odd admixture in the 
couplings of the dominantly CP-even non-SM-like light Higgs boson may weaken its couplings to 
the weak gauge bosons to such an extent that it escapes the LEP limits 110J. Additionally, CP 
violation in the Higgs couplings can allow for decays of Higgs bosons into Higgs and gauge bosons 
or a pair of lighter Higgs bosons in combinations that otherwise would be forbidden. These could 
then constitute additional discovery channels for some of the Higgs states. In general CP mixing 
alters the Higgs couplings and hence the production and decay rates with effects on the discovery 

1 Besides the fact, that the 125 GeV Higgs boson is too heavy to allow for the strong first order phase transition 
in the early universe required for thermal non-equilibrium. 
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prospects of the additional NMSSM Higgs bosonsj^] Moreover, CP-violating phases influence the 
Higgs mass spectrum already at tree level in case of Higgs sector CP violation and at loop level 
through radiatively induced CP violation 12,131. Some of these effects, however, can also be 


explained in the CP-conserving NMSSM by choosing the parameter combinations accordingly. 
Additionally, without further information on the CP nature of the Higgs bosons from other 
observables, genuine CP-violating Higgs decays to gauge boson plus Higgs or Higgs-to-Higgs 
decays cannot be identified unless all of these decay channels are observed. Since the absolute 
size of CP violation is stringently constrained by experiment, the identification of CP-violating 
Higgs bosons will be a non-trivial task requiring precision measurements, high luminosities and 
the combination of various CP-violating observables. In particular tight constraints on the CP- 
violating phases arise from the non-observation of electric dipole moments (EDMs). Accordingly, 
we discuss in detail in the main part of the paper the role of EDMs in constraining the parameter 
space of the CP-violating NMSSM. By taking into account the latest experimental constraints 
from the Higgs data and the measurements of the EDMs, we investigate the size of the EDMs 
as a function of the CP-violating phases. This allows for conclusions on the overall allowed size 
of CP violation in the NMSSM in view of the newest experimental results. Subsequently, we 
investigate how CP violation can be identified by combining Higgs-to-Higgs decays and Higgs 
decays to gauge boson plus Higgs. This is complemented by the discussion of measuring CP 
violation in fermionic Higgs decays. 

The investigation of the CP-violating NMSSM Higgs sector considered here takes into account 
higher order corrections both to the parameters and the observables. Radiative corrections to 
the Higgs boson masses are crucial to lift the SM-like Higgs boson mass to the observed value 
of 125 GeV. The connection of the Higgs self-couplings and masses through the Higgs potential 
requires the inclusion of higher order corrections also in the Higgs self-interactions to consistently 
describe Higgs-to-Higgs decays, which can become relevant for spectra with light Higgs bosons 
and/or sizeable values of the singlet coupling A [l4-24j. In the Higgs boson decays into SM 
particles and coloured SUSY particles in particular QCD corrections play an important role 
and have to be included. Also electroweak corrections can become important. They cannot be 
adapted from the SM or MSSM, however, but require an explicit calculation, as has been done 
so far in the NMSSM only for the Higgs boson decays into squarks |25]. 

It is worth recalling the status of higher order corrections to the NMSSM Higgs boson masses. 
In the CP-conserving NMSSM the one-loop mass corrections are available |26-jj35], as well as 
two-loop results of 0(cq<a s + cq><a s ) in the approximation of zero external momentum 33 . First 
corrections beyond order O (atot 8 + oibCt s ) have been provided in 1361. The one-loop corrections 
to the Higgs masses of the CP-violating NMSSM have been calculated by PL2,37-EI and the 


logarithmically enhanced two-loop effects have been given in [42]. The two-loop corrections to the 
Higgs boson masses of the CP-violating NMSSM in the Feynman diagrammatic approach with 
vanishing external momentum at 0(cq<a s ) have been computed in |13j|. The one-loop corrections 
to the trilinear Higgs self-couplings for the CP-conserving NMSSM |15j are available and have 
recently been extended to include the two-loop corrections at order 0(ata s ) in the approximation 
of vanishing external momentum in the CP-violating case |43|. Several public codes are on the 
market for the computation of the NMSSM mass spectrum for a given parameter set. These are 

and NMSSMCALC [49,50], and the 
52] and SPheno 1 53, 


the stand-alone codes NMSSMTools 144-461, SOFTSUSY 47 


NMSSM implementations of the programs FlexibleSUSY 


48 

51 


2 For a recent discussion on two Higgs bosons near 125 GeV in the complex NMSSM, see 
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based on SARAH (36,55-58] J^] An extension of the program package NMSSMTools to include also 
the CP-violating NMSSM has recently been announced in 1601. 


For the phenomenological analysis in this paper we have implemented the EDMs in the 
Fortran package NMSSMCALC. It computes for the CP-conserving and CP-violating case the two- 
loop NMSSM Higgs boson masses at 0(a t a s ) and the Higgs boson widths and branching ratios 
including the dominant higher order corrections^] 


The layout of the paper is as follows. In Section [2] we introduce the CP-violating NMSSM 
Lagrangian and set our notation. Section [3] briefly recapitulates the computation of the SUSY 
contributions to the EDMs. In Section [4] we discuss in detail the EDMs induced by various CP- 
violating phases present in the NMSSM. In subsection [4.1| this investigation is performed in the 
subspace of the Natural NMSSM that features a rather light Higgs spectrum with good discovery 
prospects for all Higgs bosons of the NMSSM. In subsection |4.2[ we extend the investigation 
to an enlarged NMSSM parameter range. Section [5] is devoted to the prospects for measuring 
CP violation in the Higgs couplings involving Z bosons through the combination of signal rates, 
subsection 5T, and in the Higgs couplings to fermions, subsection Y2 Section [6] summarizes 


and concludes the paper. 


2 The Lagrangian of the CP-violating NMSSM 

We work in the framework of the CP-violating NMSSM with a scale-invariant superpotential 
and a discrete Z 3 4 symmetry. The Higgs potential is obtained from the superpotential, the soft 
SUSY breaking Lagrangian and the D-term contributions. The NMSSM superpotential in terms 
of the two Higgs doublet superfields H^ and H u , the singlet superfield S', the quark and lepton 
superfields and their charged conjugates, with the superscript c, Q , U c , D c , L, E c , is given by 

Wnmssm = eij[y e H d L 3 E c + y d H l d Q j D c - y u H l u Q j U c ] ~ eijXSH l d H 3 u + L,S' 3 . (2.1) 

The z, j = 1,2 are the indices of the SU(2)l fundamental representation, and is the totally 
antisymmetric tensor with ei 2 = e 12 = 1, where we adopt the convention to sum over equal 
indices. Colour and generation indices have been suppressed. As we neglect generation mix¬ 
ing, the Yukawa couplings y e ,yd and y u are diagonal, and complex phases can be reabsorbed 
by redefining the quark fields without effect on the physical meaning |9|. The dimensionless 
parameters A and k, are complex in case of CP violation. 

The soft SUSY breaking Lagrangian in terms of the scalar component fields H U ^H^ and S 
reads 

£soft, NMSSM = - rn 2 H d H d H d - m 2 Hu HlH u - rn^Q^Q - m\Di - m\ R u R u R - m\Jfr R d R 

- ml R e* R e R - (e^A^L 3 S* R + y d A d H l d Q 3 d* R - y u A u H l u Q 3 u* R ] + h.c.) 

- + M 2 WiWi + M 3 GG + h.c.) (2.2) 

- m 2 s \S\ 2 + (e t] \A x SW d Hi - ^kA k S 3 + h.c.) , 

3 For a comparison of the codes, see J59 . 

4 The program package NMSSMCALC including the computation of the EDMs is made publicly available and can 
be downloaded from the url: http://www.itp.kit.edu/~maggie/NMSSMCALC/ 
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where a sum over all three quark and lepton generations is implicit. The Q and L denote the 
complex scalar components of the corresponding quark and lepton superfields and, e.g. for the 
first generation are Q = and L = In the CP-violating NMSSM the soft 

SUSY breaking trilinear couplings A x (x — A,ft,d, u,e) and the gaugino mass parameters M 
(k = 1, 2, 3) of the bino, wino and gluino fields 5, W{ (i = 1, 2, 3) and (5, are complex. Exploiting 
the R-symmetry either M\ or M 2 can be chosen to be real. On the other hand the soft SUSY 
breaking mass parameters of the scalar fields, m 2 x (X — S', H u , Q, ur , cIr , L, e#) are real. 

The Higgs potential finally is obtained as 

V H = (\\S\ 2 + m% d )HlH d + (\\S\ 2 + mljHtH u + m 2 s \S\ 2 

+\{g 2 2 + gl){H\H d - HlH u ) 2 + \gl\H\H u \ 2 (2.3) 

+ | - jj\H dti H uJ + kS 2 I 2 + [ - e lJ XAxSH d l H UJ + I k A k S 3 + h.c.] , 

with gi and g 2 denoting the U(l)y and SU(2)r gauge couplings, respectively. Expanding the 
two Higgs doublets and the singlet field about their VEVs, v d, v u and v s , two more CP-violating 
phase, (f u and tp s , are introduced, 


Hr 


_ (T hd ~\~ ICLd) 

= l h- 


H v = e Z( ^ 



+ 


S 


e i^s 

— i=( v s h s + ^a s ) • (2-4) 

v2 


The VEVs and are related to v ~ 246 GeV through v 2 = v\ + and their ratio is 
parametrized by tan/3 = v u /vd . The phase affects the top quark mass. We absorb this phase 
into the left-handed and right-handed top fields through 

tL ~^ 6 ^ u ^ 2 tx, and tR e^ u ^ 2 tR , (2-5) 


so that the top Yukawa coupling is kept real. This alters all couplings with one top quark. 


Replacing Eq. (2.4) in Eq. (2.3) yields the Higgs potential 


Vh — V£ onst + th H hd + th u h u + th s h s + taidd + t au a u + t as a 


L a s ^s 


( 2 . 6 ) 


+ -(/) O,T Al00 (f)° + c/f^Mh+h- ft 2 + Vr ^ j 


with 0° = (/i^,/i u ,/i s , a^, a n , a s ) T and <p c = ((/i rf )*,/i+) T . The tadpole coefficients are denoted 
by (</> = /^, h u , /i s , a^, a n , a s ), is the 6x6 mass matrix for the neutral Higgs bosons and 

Mh+h~ fhe 2 x 2 mass matrix for the charged Higgs states. Constant terms are summarized in 
yg° nst an( j trilinear and quartic interactions in V^ . A few remarks on the tadpoles and 
mass matrices are in order, without repeating their explicit expressions here, which can be found 
in Ref. 
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The tadpole coefficients vanish at tree level due to the minimization conditions of 
the Higgs potential. Rewriting the complex parameters A, ft, A\ and A K as 


A = \X\e l(px , ft = |ft|e zc ^ , A\ = \A\\e tipA * and A K = \A^\e lifA,<i , (2.7) 

three phase combinations appear at tree level in the tadpoles and the mass matrices, 

g>x = + <P1, (2.8) 

Vy = <^2 - <fi , (2.9) 

Vz = Va k + ¥>2 , ( 2 . 10 ) 
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where we have introduced 


+ Ps + Pu 5 (2.11) 

(f2 = <Pk H - . (2.12) 


Two of the three combinations Eqs. (2.8)-(2.10) can be eliminated at lowest order by applying 
the minimization conditions t ad = t as — 0. We express cp x and (p z in terms of cp y . All mass matrix 
elements mixing the CP-even and CP-odd interaction states, Mh iaj i are then proportional to 
sin^y. At tree level, this is the only CP-violating phase in the Higgs sector. The rotation from 
the interaction to the mass eigenstates hi (i = 1,..., 5) is performed by applying two consecutive 
rotations. The first rotation with matrix 1Z G separates the would-be Goldstone bosons. The 
second one with the matrix 1Z performs the rotation to the mass eigenstates, i.e. 


(hdi hui hs . a, a s , G) 
(hi,h 2 ,h 3 ,h4,h 5 ,G) T 

— 7 Z (J^di h U j h s , cl u , o, s ) , 

= U (h d ,h u ,h s ,a,a s ,G) T , 

(2.13) 

with the diagonal mass matrix 



diag (m 2 hl , m| 2 , m 2 h3 , m 2 hA , m 2 hh , 0) 

= UMhh'R' 7 ' , M-hh = M-4>4>('R G ) T ■ 

(2.14) 


The mass eigenstates hi are ordered by ascending mass, where the lightest mass is given by 


The tree-level Higgs potential can be parametrized by the following set of independent pa¬ 
rameters 

(2.15) 


th d , th u , th s , ta d , ta s , M 2 h± , v , sin 9 W , e, tan /3, | A |, v s , | k |, Re A K , sin ip y . 


We have chosen to use v and sin 0w, where 0w denotes the weak mixing angle, instead of Mw 
and Mz- This is more convenient in view of the inclusion of the two-loop corrections to the 
Higgs boson masses in the gaugeless limit 0 Furthermore, in accordance with the SUSY Les 
Houches Accord (SLHA) |61,62 the real part of A K is used as input parameter. The imaginary 
part is eliminated by the tadpole conditions. This distinction is not necessary for A and ft, as 
both the real and imaginary parts are given in the SLHA convention and can be related to the 
respective absolute values and phases. Note finally, that the effective higgsino mixing parameter 
is given by 


Meff = 


|A| y s e 


i&s+y a ) 


V2 


(2.16) 


At higher order in the Higgs mass corrections (and self-couplings) the CP-violating phases 
entering the Higgs sector at tree level are not related any more. The phases cpi and cp 2 appear 
independently in the neutralino sector, while the chargino and up-type squark sector depend on 
the phase Lp\. We therefore have two independent CP-violating phases that appear, if we choose 
to determine the phases (pA K and cpA x from the tadpole conditions. Of course non-vanishing (p\ 
or will automatically imply non-vanishing ^PA\ and However, since we consider the latter 
two as derived quantities, which are fixed via the tadpole conditions, we will not mention them 
explicitly in the following discussion. The higher order corrections introduce further complex 
phases stemming from the complex soft SUSY breaking trilinear couplings and gaugino mass 
parameters, that enter the couplings and SUSY particle masses involved in the loop corrections. 

5 See |l3j for details. 
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Figure 1: Generic one-loop diagrams contributing to the EDMs of the electron and light quarks (/ = e,u,d,s). 





Figure 2: Generic one-loop diagrams contributing to the CEDMs of the light quarks (q = u , d , s). 


3 Electric Dipole Moments 

CP violation would manifest itself in the generation of EDMs. The non-observation of any 
EDMs so far poses stringent constraints on CP-violating phases. These have to be taken into 
account, when discussing possible CP-violating effects in the Higgs sector. The constraints may 
become weaker in case of accidental cancellations of the various contributions to the EDMs, 
even for lighter SUSY particles with masses below 0(1 TeV) |63|[64]. In our analysis we take 
into account all relevant CP-induced contributions to the observable EDMs. In particular we 
consider the compatibility with the experimental upper bounds on the EDMs, which are 


Electron EDM |65j 

~ 1 • 10 -28 ecm 

Thallium EDM 66 

~ 9 • 10 _25 e cm 

Neutron EDM 67 

~ 3 • 10 -26 e cm 

Mercury EDM 68 

~ 3.1 • 10 -29 ecm 


where the electron EDM is estimated from the thorium monoxide experiment. These observable 
EDMs receive form factor contributions from the electric dipole moment, the chromo-electric 
dipole moment (CEDM), the two-loop Weinberg three-gluon operator and the Higgs-exchange 
four-fermion operators. All of these contain contributions that are generated by CP-violating 
Higgs mixing at tree level. In the EDM and the CEDM we consider one- and two-loop contri¬ 
butions. The two-loop contributions stemming from CP violation in the Higgs sector at tree 
level, which is specific to the NMSSM, can become important when the CP phases of the MSSM 
parameters are set to zero 1 69 [. Such a configuration can be achieved by choosing cp2 ^ 0, while 
keeping cpi = 0. Therefore we will refer to ip 2 as the NMSSM-specific phase in the following. 
The one-loop EDMs of the electron and the light quarks u , d, s are induced by chargino, xf 
(j = 1,2), and neutralino, Xi — 1, 5), exchange diagrams, cf. Fig. [l] For the quarks also 

gluino, g, exchange diagrams contribute. The light quarks furthermore have CEDMs which are 
also generated by chargino, neutralino and gluino loops as shown in Fig. [2j The one-loop EDMs 
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Figure 3: Generic two-loop Barr-Zee type diagrams contributing to the EDMs of the electron and light quarks 
(/ = e,u,d,s). Upper: 7 i7, and lower from left to right: ZH , WH and WW. 



have been computed before and the formulae can be found in [64j[69j. At two-loop level the 
Higgs mediated Barr-Zee type diagrams contribute significantly to the EDMs. They are medi¬ 
ated by neutral Higgs couplings to two photons, 77 Hf [64j, the charged Higgs coupling to the 
charged W boson and a photon, r yH ± W T , the r yW ± W T coupling [70|, and the couplings between 
a neutral Higgs boson, a photon and a Z boson, 7 HfZ |711. We will denote these contributions 
in the following as yi7, WH , WW and ZH , respectively. The diagrams are displayed in Fig. [3) 
Additionally, CEDMs of the light quarks iq d and s are generated by two-loop Higgs-mediated 
Barr-Zee graphs 1 64,69 , cf. Fig. [4j For the Weinberg operator we take into account the con¬ 


tributions from the Higgs-mediated two-loop diagrams \72\ and additionally the contribution 
from the quark-squark-gluino exchange contribution 173]. The coefficients of the four-fermion 
operators, finally, are generated from the ^-channel exchanges of the CP-violating neutral Higgs 
bosons [64] . 

We briefly describe how the observable EDMs are obtained from the various contributions 
introduced above. For explicit formulae and more details we refer the reader to |64,69| and 
references therein. The Thallium EDM receives contributions from the electron EDM and ad¬ 
ditionally from the CP-odd electron-nucleon interaction [74, 75 . For the neutron EDM three 


different hadronic approaches are considered. These are the Chiral Quark Model (CQM), the 
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Parton Quark Model (PQM) and the QCD sum rule technique. In the non-relativistic CQM 


the quark EDMs are estimated via naive dimensional analysis 76 . The PQM uses low-energy 


data related to the constituent-quark contribution to the proton spin combined with isospin 
In the third approach finally QCD sum rule techniques (78-801 are applied to 


77 


symmetry 

determine the EDM. Depending on the approach, the neutron EDM is composed of the contri¬ 
butions from the EDMs and CEDMs of the light quarks, the Weinberg operator and the CP-odd 
four-fermion operators, see (64j for details. By using QCD sum rules (78|[80| the Mercury EDM 
is estimated from the EDMs induced by the Schiff moment^J the electron EDM, the contribu¬ 
tion due to the CP-odd electron-nucleon interaction and the contributions from the couplings 
of electron-nucleon interactions. Details are given in [j69|. 


4 EDM constraints 

In this section we investigate the influence of the various CP-violating phases on the EDMs 
and the resulting constraints on possible CP-violating scenarios. We will present results for the 
phases pi,p2 and the phase pA t i which arises from the top/stop sector, and the phases pMi of 
the gaugino mass parameters Mi (z = 1,2,3). The EDMs also depend on the phases pa 3 of the 
soft SUSY breaking trilinear couplings Aj (j = 6, iq c, d, s). The influence of the phases pa- of 
the sbottom and the first and second generation squark sector is by far subleading compared to 
the effects of the other phases. The reason is that the trilinear couplings and hence their phases 
come in combination with the quark masses, which are small in this case. 

In order to find viable parameter points we perform a scan in the NMSSM parameter space 
and keep only those points that are in accordance with the LHC Higgs data. That this is the case 
has been checked with the help of the programs HiggsBounds [82}|84] and HiggsSignals |85|. 
For the computation of the Higgs boson masses, the effective couplings, the decay widths and 
branching ratios of the SM and NMSSM Higgs bosons, that are needed as inputs for HiggsBounds 
and HiggsSignals, the Fortran code NMSSMCALC [49j, 50| has been used. Besides the masses at 
two-loop level, it provides the SM and NMSSM decay widths and branching ratios including 
the state-of-the-art higher order corrections. We demanded that the valid scenarios feature a 
Higgs boson with mass around 125 GeV. With HiggsBounds we have checked whether or not the 
Higgs spectrum is excluded at the 95% confidence level (CL) with respect to the LEP, Tevatron 
and LHC measurements. The package HiggsSignals tests for the compatibility of the SM-like 
Higgs boson with the Higgs observation data. We required the p- value, that is given out, to be 
at least 0.05, corresponding to a non-exclusion at 95% CL. 


4.1 Natural NMSSM 


We performed a parameter scan in the subspace of the Natural NMSSM as defined in 21 . It 
features a rather light overall Higgs mass spectrum and gives good discovery prospects for all 
NMSSM Higgs scalars. It is characterized by 


0.6 < |A| < 0.7, \k\ < 0.3, 1.5 < tan/3 < 2.5, 100 GeV < |/i eff | < 185 GeV. 


(4.18) 


The small ft values lead to an approximate Peccei-Quinn symmetry. Note, that in this parameter 
region the second lightest of the Higgs bosons that are dominantly CP-even, is SM- litofl The 
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For the Schiff moment several approximations exist 


81 


7 The Higgs boson with mass close to 125 GeV is forced to be mostly CP-even due to the requirement to be 
compatible with the LHC Higgs data. 








soft SUSY breaking trilinear NMSSM couplings are varied in the interval 

\A X \ < 2 TeV and \A K \ < 2 TeV . (4.19) 

The remaining soft SUSY breaking trilinear couplings and masses have been chosen as 

\Au\, \Ap\, \ Al\ < 2 TeV with U = u,c,t, ,D = d,s,b,L = e,n,r, (4.20) 


Mu R ,c R — Md R ~ R — Mq i 2 — Me R ^ R — M^ i 2 — 3 TeV , 


(4.21) 


600 GeV < M 1r = Mq 3 < 3 TeV , 600 GeV < M? R = M u < 3 TeV , M Sp = 3 TeV , (4.22) 


l 3 


bR 


100 GeV < \Mi\ < 1 TeV , 200 GeV < \M 2 \ < 1 TeV , 1.3 TeV < |M 3 | < 3 TeV . (4.23) 

All scenarios have been checked for compatibility with the lower bound on the charged Higgs 

Note also, that 


86 and respect the exclusion limits on the SUSY particle masses 187-89 


mass 

the signs of the gaugino masses M\ and M 2 have only a marginal effect on the features of the 
NMSSM Higgs sector. The NMSSM-specific input parameters A, ft, A\ and A K as well as all 
other soft SUSY breaking masses and trilinear couplings, according to the SLHA format, are 
understood as DR parameters taken at the SUSY scale Ms = ^In NMSSMCALC also 
tan/3 is assumed to be given at the SUSY scale. 

Variation of (p 2 : For the scenarios of our scan that are compatible with all above described 
constraints we computed the various EDMs and checked for their compatibility with the exper¬ 
imental values. The results are shown in Fig. [5] for the electron, the Thallium, the neutron and 
the Mercury EDM. In these plots we have only varied the NMSSM specific phase <p 2 . All other 
complex phases have been set to zero. We hence investigate here solely tree-level CP violation, 
as it can occur in the NMSSM. As the phase <p 2 only appears in the NMSSM we call this in 
the following NMSSM-type CP violation. For each EDM, the absolute values of the computed 
EDMs in the valid scenarios are shown, normalized to the respective experimental upper bound, 
as given in Eq. (3.17). All points above 1 are hence in conflict with the EDM data. For the 


neutron EDM results are given for two different hadronic approaches, the chiral quark model and 
the one applying QCD sum rule techniques, Fig. [H] (middle row) J^] In the Mercury EDM, Fig. [ 5 ] 
(lower row), we have taken into account the uncertainties in the calculation of the contribution 
from the Schiff moment by showing results for two of the four values given in the literature, 

dip, m , IV [s] 


;n g 

inferred from t 




69,8l|. The superscripts in da g indicate which value has been applied. As can be 
ie"plots the most stringent constraints arise from the electron EDM, Fig. [ 5 ] (upper 
left). The constraints from the Thallium EDM, Fig. [5] (upper right), which mostly depends on 
the electron EDM dE e , are not as stringent. The results for the neutron EDM based on the 


e technique differ by about a factor 2-3 in accordance with previous 


CQM and the QCD sum ru 
results in the literature 60,69 . The pictures for the Mercury EDM based on d^ g and d^ are 

almost the same and similar to the ones based on d^ and d^J, which are not shown here. The 
electron EDM is larger for larger values of |ftA|, i.e. |ftA| > 0.1. Note, that CP violation in the 
tree-level Higgs sector is induced by terms proportional to |ftA|sin(<^i — cp 2 ) so that for larger 
absolute values of |ftA| the CP-violating phase plays a more important role in the EDM. The 
figures show, that despite the exclusion of some points by the tight electron EDM constraints, 
scenarios are viable that feature a large CP-violating phase, including possible maximum CP 
violation cp 2 — ±7t/ 2. Finally let us remark that the asymmetry in the scattering of the points 
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Figure 5: Absolute values of the electron (upper left), Thallium (upper right), neutron (middle) and Mercury 
(lower) EDMs as a function of (^ 2 , normalized to the respective experimental upper bound; red: \k,\\ < 0.1, blue: 
\kX\ > 0.1. 


at cp 2 = 0 and ±7r is due to the fact, that our scan only extends over positive values of A and 
that k, appears always in the product ftA. 

Investigating the various contributions in detail, we find that at one-loop level only the 
neutralino exchange diagrams contribute to the electron EDM. All other one-loop contributions 
do not depend on cp 2 . The one-loop contributions turn out to be well below the exclusion bounds. 
The dominant part comes from the two-loop diagrams where the jH contribution is the most 
relevant one, but also WW is significant and comes with a different sign, see Fig. [6j The WH 
and ZH contributions are about an order of magnitude smaller. 

The analysis of the Thallium EDM shows that the one-loop contributions are tiny. The 
8 In NMSSMCALC we also implemented the Parton Quark Model, but do not show explicit results here. 
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two-loop yH 
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Figure 6 : Absolute values of the two-loop 7 H (left) and WW (right) contributions to the electron EDM normalized 
to the experimental upper bound. Red points represent contributions with negative sign, green points those with 
positive sign. 


two-loop and the four-fermion operator contributions are roughly of the same size and come 
with opposite sign, so that they partly cancel each other. Concerning the neutron EDM both 
in the CQM and in the QCD sum rule approximation the one-loop contributions are irrelevant 
and the two-loop part is dominated by the Weinberg operator contribution and the Barr-Zee 
type contributions to the CEDM of the quarks. In the QCD sum rule approach the latter two 
come with opposite signs. The four-fermion operator part that also contributes here, is small. 
The largest contributions to the Mercury EDM originate from the electron EDM. Also the 
contributions induced by the Schiff moment are of comparable order, but come with a different 
sign. 

Variation of ipA t •' We now turn to the discussion of the effects of a non-vanishing phase of 
the stop trilinear soft SUSY breaking coupling, (pA t - This is a phase that also appears in the 
MSSM and we call this in the following MSSM-type CP violation. Figure [7] shows the absolute 
values of the observable EDMs computed for the NMSSM scenarios from our scan normalized 
to the experimental upper bounds, now as a function of (fA t - All other possible CP-violating 
phases have been set to zero. As can be inferred from the figure, the most important observable 
EDMs induced by ipA t are the electron EDM dE e , Fig. [T| (upper left), and the neutron EDM 
dn, Fig. [7] (middle row), both, however, being well below the experimental values for most of 
the parameter points. In the Mercury EDMs differences in the application of different Schiff 
moment contributions are now well visible. They are not of relevance though, as the obtained 
values remain below the experimental limit. 

Note that a non-zero phase ipA t can induce at one-loop level only for the up-quark an EDM 
and a CEDM through the neutralino and gluino exchange diagrams, but solely when stops are 
involved, cf. Figs, [l] and [ 2 | At two-loop level, besides the Weinberg operator, only the Barr-Zee 
type diagrams jH lead to a non-vanishing contribution, cf. Figs. [3] and [4j All other diagrams 
involve no couplings which contain ipA t . Therefore the electron and the Thallium EDM only 
receive two-loop 7 # contributions, cf. Fig. [3] (upper), which change sign with the phase (pA t - 
In the neutron EDM the two-loop contributions, which also include those from the Weinberg 
operator, are about one to two orders of magnitude larger than the one-loop contributions, 
depending on the approximation. Both come with opposite sign, cf. Fig. | 8 j The Weinberg 
operator provides the dominant part at two-loop level, while the EDM and CEDM contributions 
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Figure 7: Same as 


Fig. 1 but for the variation of ipA t • 


All other CP-violating phases are set to zero. 


to the quarks, which enter at one- and two-loop level, are of about the same size and come with 
opposite sign. 


Variation of <pi: A non-vanishing phase pi leads to CP violation in the Higgs sector already 
at tree level and hence generates an NMSSM-type CP violation. As the phase <pi also enters 
the effective higgsino parameter fi e fp cf. Eq. (2.16), it also generates CP violation in the doublet 
higgsino and in the sfermion sector as it occurs in the MSSM, and therefore leads to MSSM-type 
CP violation. In Fig. [9] we show the effect of a complex phase <pi on the values of the EDMs. 
All other phases have been set to zero. As can be inferred from the plots, CP violation induced 
by cpi is strongly constrained by the EDMs. In particular the induced electron EDM is by a 
factor up to 100 times larger than the experimental upper bound. At one-loop level the chargino 
contributions, that are MSSM-like, are important. The Barr-Zee type two-loop contributions 
are an order of magnitude larger than the ones induced by cp 2 ^ 0. In the Thallium EDM the 
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Figure 8: Absolute values of the one- (left) and two-loop contributions to the neutron EDM in the QCD sum 
rule approach, normalized to the measured upper bound. Red points represent contributions with negative sign, 
green points those with positive sign. 



Figure 9: Absolute values of the electron (upper left), Thallium (upper right), neutron (middle) and Mercury 
(lower) EDMs as a function of <^i, normalized to the respective experimental upper bound. 
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one- and two-loop contributions are of same size and come with the same sign. They dominate 
over the four-fermion operator part. For the neutron EDM the one-loop contributions dominate, 
while the contributions originating from the two-loop diagrams and the Weinberg operator are 
about one order of magnitude smaller and cancel each other partly. In the Mercury EDM, the 
dominant part is built up by the electron EDM. 

In order to disentangle how much of the observed effect originates from NMSSM-type or, 
respectively, MSSM-type CP violation, we varied the phases ip\ and at the same time. Setting 
p i = ip 2 allows to turn off tree-level CP violation in the NMSSM Higgs sector, since this yields a 
vanishing CP-violating phase in the tree-level Higgs sector. Hence only the CP-violating MSSM- 
type contributions induced by p\ remain. The resulting plots, which we do not show explicitly 
here, look strikingly similar to Fig. [9} The major difference is that in Fig. [9]there are hardly any 
points around the maximally CP-violating phases cpi — ±7t/ 2, whereas for the cpi = cp 2 variation 
these points exist. The gap in Fig. [9] around p\ — ±7r/2 can be attributed to the fact that we 
demand compatibility with the Higgs data. However, if the CP violation in the tree-level Higgs 
sector is too strong, i.e. if the CP-odd admixture of the SM-like boson becomes too large, the 
signals of the SM-like Higgs boson are not compatible with the experimental values any more. 
Regarding the EDMs we conclude that in the parameter space of the Natural NMSSM the 
MSSM-type CP-violating contributions induced by p\ dominate by far over the NMSSM-ones 
generated by p\. 



Figure 10: Absolute values of the neutron (upper) and Mercury (lower) EDMs as a function of Pm 3 , normalized 
to the respective experimental upper bound. 


Variation of the phases PmS^ — 1,2,3,); We comment here on the influence of the phases 
if Mi of the gaugino mass parameters Mi on the EDMs. They are other examples for phases, 
that appear in the MSSM, too. Compared to the results from the previously discussed phases 
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(^i, if 2 and ipA t: the gaugino phases do not lead to more stringent constraints nor add new effects 
on the EDMs. 


The phase ifM 3 , that arises in the stop and gluino sector, does not contribute to the electron 
EDM. Its contributions to the neutron and Mercury EDM are of the same size as those arising 
from ifA t i as can be inferred from Fig. 10, It shows the absolute values of the neutron and the 
Mercury EDMs as a function of (fM 3 , normalized to the respective experimental upper bounds. 
Here and in the following plots, for the neutron EDM results are given in the chiral quark 
model approach and in the one based on QCD sum rule techniques. For the Mercury EDM the 
presented results are obtained for two different values of the Schiff moment, and 

as defined in 169,81. 


The phase f>M 2 plays a role in the chargino sector, just as the phase of the effective /i 



Figure 11: Absolute values of the electron (upper left), Thallium (upper right), neutron (middle) and Mercury 
(lower) EDMs as a function of (fM 2 , normalized to the respective experimental upper bound. 
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parameter, encoded in tpi. The sizes of the EDMs due to <pm 2 can therefore be expected to be 
of the same order as those due to This is confirmed by the plots given in Fig. HU which 
show the absolute values of the electron, Thallium, neutron and Mercury EDMs as a function 
of normalized to the respective experimental upper bounds. Since the phase (pM 2 has only 
a marginal effect on the Higgs sector, the compatibility of the Higgs data is not affected by (£m 2 
and no gap around ±7r/2 arises as in the results for non-zero <^i, cf. Fig. [9j 

Finally, the phase (fM 1 arises in the neutralino sector and therefore only gives contributions 
to diagrams that involve neutralinos in the loops. The generated EDMs are smaller than the 
ones for a non-vanishing (fM 2 , cf. Fig. [l2j The phase ipM 1 therefore does not lead to more 
stringent constraints on the CP-violating NMSSM as the ones that have already been discussed 
above. 



Figure 12: Absolute values of the electron (upper left), Thallium (upper right), neutron (middle) and Mercury 
(lower) EDMs as a function of normalized to the respective experimental upper bound. 
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4.2 Enlarged NMSSM Parameter Space 

We now turn to the discussion of the EDMs in an enlarged NMSSM parameter space. While 
the previous scan focused on NMSSM regions that lead to an overall light Higgs mass spectrum 
with good discovery prospects for all Higgs bosons, we here cover a large part of the NMSSM 
parameter space. In particular we also allow now for large values of tan (3 and cover the whole 
allowed space of A and n while taking care of the perturbativity constraint 

y/\X\ 2 + |k| 2 < 0.7 . (4.24) 

Additionally, we vary the effective fi parameter in a large range. In summary, our scan covers 


1 < tan/3 < 30 , |A| < 0.7 , 


\n\ < 0.7 


|/Aff| < 1 Te V 


(4.25) 


The ranges of the remaining parameters are the same as given in Eqs. (4.19)-(4.23). And we 
again checked for the compatibility with the lower bound on the charged Higgs mass 86 and 


the exclusion limits on the SUSY particle masses |86f|89| . We now also have scenarios where the 
lightest of the mostly CP-even-like Higgs bosons can be SM-like. 

Variation of (p 2 : For the scenarios of the large scan we show the results for the EDMs com¬ 


pared to the experimental values in Fig. [T3j They are plotted for the electron, the Thallium, the 
neutron and the Mercury EDM. We have only varied the NMSSM specific phase and set all 
other complex phases to zero. The comparison with the results of the Natural NMSSM in Fig. [5] 
shows, that more parameter sets now lead to EDMs that exceed the experimental limits. This is 
particularly striking for the neutron EDM, where the EDMs can be up to a factor of 20 larger. 
At the same time we have more EDMs with very small values, in particular in the case of large 
CP-violating phases. The orange points indicate the results for tan/3 values below 5, while red 
ones refer to larger tan /3 values. This shows the strong influence of the tan /3 parameter on the 
size of the EDMs. Larger tan /3 values lead to larger EDMs and vice versa. This is in accordance 
with the findings of Refs. 64,90]. Note, finally that Fig. 13 and all subsequent figures of this 


subsection of course contain the subspace of the Natural NMSSM. As we scan here, however, 
over a larger parameter space the distribution of the points may not be exactly the same as in 
the corresponding figures for the Natural NMSSM. 

In the enlarged scan, not only the electron EDM, but now also the neutron and Mercury 
EDMs can lead to stringent constraints on the parameter space. Otherwise, the investigation of 
the individual contributions shows the same pattern as in the Natural NMSSM subspace. Only 
the sign of the various contributions is not related to the sign of the phases as clearly any more 
as e.g. in Fig. [6| 

Variation of PA t : The effect of a non-vanishing MSSM-like phase pA t is shown in Fig. 


14 


Apart from the neutron EDM based on the CQM approach all EDMs display enhancements, 
sometimes by up to a factor of 50 for the maximal values, due to the large parameter space. 
This is in particular the case for small phase values, \pA t \ ^ 0.25 7r, and large values of tan/3, 
where now all EDMs contribute to the exclusion bounds. 

The individual features of the EDMs as discussed already for the Natural NMSSM, do not 
change. 

Variation of : Finally, we show the constraints arising from a non-vanishing phase p\. The 
enlargement of the parameter space leads to an increase of the maximal values of all EDMs by 
about a factor of 10, cf. Fig. [15) The investigation of the EDMs shows, that for the variation 


of pi they are very sensitive to the value of tan/3. This can be inferred from Fig. 15 where the 
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Figure 13: Large parameter scan: Absolute values of the electron (upper left), Thallium (upper right), neutron 
(middle) and Mercury (lower) EDMs as a function of (^ 2 , normalized to the respective experimental upper bound; 
orange: tan [3 < 5, red: tan/3 > 5. 


results for tan/3 < 5 are indicated by the orange and those for tan/3 > 5 by the red colour. 
Overall, the non-vanishing phase cpi leads to the strongest constraints on EDMs. This is mainly 
due to the MSSM-specific CP violation from a complex /i parameter. 

We also checked explicitly the behaviour of the EDMs with increasing SUSY particle masses. 
As expected we observe a decoupling behaviour, i.e. the EDMs decrease with increasing SUSY 
particle masses in the loops. 
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Figure 14: Same as Fig. 13 but for the variation of (fA t - All other CP-violating phases are set to zero. 


5 Phenomenological Higgs Analysis 

The explicit verification of CP violation at the LHC is a non-trivial task and requires high 
luminosities 91 . For the measurement of CP violation, observables can be constructed that are 
sensitive to CP-violating effects in the Higgs couplings to gauge bosons and fermionsj^] While the 
Higgs couplings to massive gauge bosons project on the CP-even component of the Higgs state, 
the fermionic couplings have the advantage to democratically couple to the CP-even and CP-odd 
components of the Higgs bosons. The individual signal rates, on the other hand, do not allow 
for conclusions on CP violation, since specific parameter configurations in the CP-conserving 
NMSSM can lead to the same rates as in the CP-violating case within the experimental errors. 


y For a comprehensive list of the relevant literature, see e.g. 91,92 and references therein, complemented by 


recent investigations in 93 113 
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Figure 15: Same as Fig. 


13 


but for the variation of ip±. All other CP-violating phases are set to zero. 


Also the superposition of rates stemming from CP-odd and CP-even Higgs bosons, that are close 
in mass, can mimic CP-violating effects. The simultaneous measurement of Higgs decay rates, 
however, that require a CP-even, respectively, a CP-odd component, allow for conclusions on 
CP violation. 

In the following two subsections we will discuss possible prospects for accessing CP violation 
in two different approaches. These are on the one hand the combined measurement of signal 
rates and on the other hand the exploitation of decays into fermion pairs. 

5.1 Hints towards CP violation in Higgs decays involving a Z boson 

One possibility to pin down CP violation in the Higgs sector is the observation of a combination 
of decays that is not allowed in CP-conserving scenarios. If a new scalar resonance is discovered 
there are several decays that offer insights on its CP nature. For example if the new resonance 
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H' decays into a pair of vector bosons, it has to have a CP-even admixture. Furthermore one 
can make use of the knowledge that the SM-like Higgs boson at 125 GeV is mostly CP-even. 
Hence, the decay H' —>► hh is also an indication for a large CP-even component of H '. However, 
if at the same time H' hZ can be observed, H r must feature a CP-odd component as well. 
In short, if we observe 


H' ZZ or H f hh = CP H / = +1 ( . 

and H' —» Zh = CP^/ = -1 

at the same time this proves CP violation in the Higgs sect or 

To illustrate this idea we present an example scenario, taken from the scan we performed 
with NMSSMCALC as described in Sec. |4j The parameters of the Higgs sector are given b)0 


|A| = 0.635 , |k| = 0.288 , \A K \ = 210.70GeV , |/x eff | = 178.65 GeV , 

<£>1=0, <£>2 = 0.0119 7T , ipA x = 0.0037 7T , (fA K = 0.983 7T , 

tan/3 = 1.88, M H ± = 392.9GeV . (5.27) 

The other input parameters are 

m u R ,CR = m d R ,s R = m Q 1>2 = m L 1:2 = m e R J>R = 3 TeV, m~ tR = 844 GeV, 
rriQ 3 — 844 GeV, rn j )R — 3 TeV, mi 3 — 1751 GeV, mf R = 1751 GeV, 

|Awl = 603 GeV, \A d:S:b \ = 16 GeV, \A e ^ r \ = 1086 GeV, (5.28) 

|Mi| = 764GeV, |M 2 | = 756GeV, |M 3 | = 2650GeV, 
a,c,< = n = VAe^r = V Mi = ( PM 2 = PM 3 = 0 • 

The resulting Higgs spectrum is relatively light, 

M Hl w 104 GeV , Mh 2 w 126.4 GeV , M Hs « 258 GeV , (5.29) 

M Ha « 401 GeV , M Hb w 405 GeV . (5.30) 


Possible candidates for H r are either H 3 , H 4 or H$. However, the two heavy mass eigenstates 
are very close in mass, so that their individual signals cannot be disentangled. In most scenarios 
of the CP-conserving case one of the heavy mass eigenstates is CP-even and the other CP-odd. 
But since it is impossible to tell whether the decay of H 4 or H$ is observed, no conclusion about 
CP violation can be drawn. This leaves us with H f = H%. For H% both the decay into a pair of 
vector bosons as well as the decay into a Z boson and the SM-like Higgs boson H 2 are sizeable, 


BR (H 3 -a ZH 2 ) = 5.7% , BR (H s -A ZZ) = 11.8% , BR (H 3 -A WW) = 27.5% . (5.31) 


These decays can only be observed if the production cross section for H% is sufficiently large. 
At y/s = 13 TeV the production cross section through gluon fusion at next-to-next-to-leading 
ordei 12 is = 211.2 fb, which is not particularly large for a Higgs boson of this mass. The 

reason is the suppression of the effective coupling of H 3 to a pair of gluons due to the large 
singlet admixture to the mass eigenstate. In Table [l] we list the obtained signal rates for several 


114 


10 For a recent investigation in the complex 2-Higgs-Doublet Model, see 

11 Note, that (fA x and (fA K are strictly speaking not input parameters, but derived quantities, determined via 
the tadpole conditions. 

12 The cross section has been calculated with a private version of HIGLU 115 , that has been adapted to the 
complex NMSSM. 
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t 


a(ggH 3 ) BR (H 3 ZZ) 

24.8 fb 

a(ggH 3 ) BR (H 3 -)■ WW) 

58.1 fb 

<r(ggH 3 ) BR (H 3 ->■ H 2 Z ) 

12.1 fb 

a{ggH 3 ) BR (H 3 ->• H 2 Z ->• (bb)Z) 

5.5 fb 

<r(ggH 3 ) BR (H 3 ->■ H 2 Z ->■ ( 77 )Z) 

0.04 fb 

a(ggH 3 ) BR (H 3 ->■ H 2 Z ->■ (ZZ)Z) 

0.45 fb 

a{ggH 3 ) BR (H 3 H 2 Z (WW)Z) 

3.5 fb 

a(ggH 3 ) BR(H 3 ->■ H ± Z) 

76.0 fb 

a(ggH 3 ) BR (H 3 ->■ H X Z ->■ (bb)Z) 

66.0 fb 

a(ggH 3 ) BR (H 3 ->• H 3 Z ->• ( 77 )Z) 

0.04 fb 

< 7 ( 55 ^ 3 ) BR(F 3 -)■ HiZ -)■ (ZZ)Z) 

0.06 fb 

£7( 55 j ff 3 ) BR(iJ 3 -> HaZ (WW)Z) 

0.65 fb 

BR(ff! ->• ZZ) 

1.53 fb 


Table 1: Signal rates for H 3 production in gluon fusion at y/s = 13 TeV with subsequent decay into various final 
states. 


final states. The rates in the vector boson final states ZZ and WW are with 24.8 fb and 58.1 fb 
of moderate size. For the final state ZH 2 the signal only amounts to 12.1 fb. If the decay of the 
SM-like H 2 is taken into account as well the overall signal becomes rather small. The decay of 
H% into a Z boson and Hi, on the other hand yields a signal rate of 76 fb. Together with the 
decay of the lightest Higgs boson into a pair of Z bosons, which establishes a CP-even admixture 
to Hi, this could also be used to search for CP violation in the Higgs sector. The signals are 
rather small though so that the experimental observation will be challenging. 

Note that in this scenario the EDM bounds are respected. As we observed in Sec. [4j this is 
not uncommon for scenarios that feature (^2 ^0 and cpi = 0. 


5.2 CP violation in the Higgs to r + r decays 


In 103,116 118 and 94 the possibility to make use of the r + r _ decay mode for the determi¬ 
nation of the CP properties of a Higgs boson has been suggested. We therefore investigate here 
the coupling of Higgs bosons to a r pair. If a mass eigenstate Hi is an admixture of CP-even and 
CP-odd components, it can feature both a scalar coupling (denoted by cf) and a pseudoscalar 
coupling (denoted by c^). The corresponding term in the Lagrangian reads 


jC-HiTT = - 9 t T + (cf + icf 75 ) T Hi 


m T 


with g T =- 


(5.32) 


f - _ c p 

— ~9r\J ( c r) 2 + ( c t ) 2 r+ ( cos Pi + i sin <^ 75 ) r Hi with tan<^ = ^r. (5.33) 


In the second line we followed 1103,119j and introduced the CP-violating angle (pi which parame¬ 
trizes the CP-mixing of the Higgs boson Hi which couples to the r. For a CP-even Higgs boson 
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Figure 16: The phase 0$, which measures the CP violation in the Hir + r~ coupling, as a function of the mass 
of the Higgs boson Hi. Only scenarios that feature H 2 as the SM-like Higgs boson with a mass of 125 GeV are 
included. Light colored open points are in conflict with the EDM bounds, whereas dark colored full points respect 
the bounds. 


4>i = 0, and for a CP-odd one <pi — 7 t/2 . Figure 16 shows pi plotted against the mass of the 
respective Higgs boson. We included all points of the previously described scans that feature the 
next-to-lightest Higgs boson as the 125 GeV Higgs state. The light colored open points indicate 
conflict with the EDM constraints, while dark colored full points respect the EDM constraints. 
Figure [16] shows that H\ and H 2 couple mostly scalar-like to the r pair, whereas H 3 features a 
pseudoscalar-like coupling in most cases. For the heavier Higgs states no such tendency can be 
observed. Especially for the three heavier Higgs states there are points that display CP violation 
in the HiT + r~ coupling and that are not in conflict with the EDM bounds. 


In the following we investigate by which phases of the input parameters these CP-violating 
effects in the HiT + r~ coupling can be generated. Figure pT| shows the phase pi plotted versus 
several phases of the input parameters, namely (upper left), ^2 (upper right), ipA t (lower left) 
and cpi = (f 2 (lower right). Even relatively small phases of ipi and cp 2 can generate considerable 
CP violation in the HiT + r~ couplings that should be accessible in the experiment. This is due 
to CP violation already at tree level in the Higgs sector, that can occur if solely or cp 2 are 
set to non-trivial values. As we already saw earlier the EDM bounds are exceeded, however, 
if (fi takes on non-trivial values, whereas c P 2 leads to scenarios, which respect the bounds. If 
the CP violation is only induced by loop effects, as it is the case for cpA t or if cpi and cp 2 are 
chosen equal, hardly any CP-violating effect is visible in the i^T + r _ coupling. The prospects 
of measuring such a small CP violation are less good in this case. 


6 Conclusions 

Supersymmetric theories feature many new sources for CP violation. In particular in the 
NMSSM, CP violation can already be induced at tree level in the Higgs sector. The upper 
bounds on the EDMs, on the other hand, pose stringent constraints on possible CP-violating 
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Figure 17: The angle 0* as a function of various phases of the input parameters. All other phases are set to zero. 
Light colored open points are in conflict with the EDM constraints, whereas dark colored full points respect the 
bounds. 


phases. We have investigated the allowed ranges for CP-violating phases in the NMSSM by 
taking into account the current limits on EDMs and the latest Higgs data from the LHC. On 
the one hand the phase v?i — ^ 2 , which induces the tree-level CP violation in the Higgs sector, 
is of interest. On the other hand radiative corrections to the Higgs masses are indispensable 
to achieve the measured mass value of 125 GeV for the SM-like Higgs boson. In particular the 
top/stop sector, which delivers the dominant Higgs mass corrections, introduces the additional 
phase (fA t - Furthermore the phase cpi appears on its own (not in combination with (^ 2 ) in the 
stop sector and also enters in the chargino sector. 

Our analysis has shown that the EDMs induced by the NMSSM specific phase (^ 2 , that 
appears already at tree level in the NMSSM Higgs sector in contrast to the MSSM, leads to small 
contributions to the EDMs. In this case, non-trivial CP-violating phases can still be compatible 
with the EDMs. The most stringent constraints on the phases are then due to the LHC Higgs 
data, as the possible large CP admixture in the 125 GeV Higgs boson leads to signal rates that 
are not compatible with the experiment any more. On the other hand chargino contributions 
to the EDMs through a complex phase of the effective fi parameter generate EDMs that are 
above the experimental constraints. The EDM contributions stemming from a non-vanishing 
phase (fA t can be important for small values of the phases and large values of tan/3. Otherwise, 
this phase is hardly constrained by the EDMs. Overall, we observe that the induced EDMs are 
more important for larger values of tan/3. Finally the phases of the gaugino mass parameters 
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can induce sizeable EDM contributions. We have also verified, that the CP-violating effects in 
the EDMs decouple with rising masses of the SUSY particles in the loops as expected. 

The experimental verification of CP violation in the Higgs sector is a non-trivial task. Higgs 
signal rates can be used for the verification of CP violation only in the simultaneous measurement 
of decay rates requiring a dominantly CP-even, respectively, CP-odd admixture in the Higgs mass 
eigenstate. We have shown one example where such measurements should be feasible, although 
the signal rates are challenging. Besides other observables, Higgs decays into fermions provide 
observables sensitive to CP violation. Taking into account the EDM and LHC Higgs constraints 
we have shown, that in particular for the three heavier Higgs bosons there are scenarios that 
display CP violation in the Higgs couplings to tau leptons and are not in conflict with the 
bounds on the EDMs. Our investigation of the origin of the CP effects has revealed, that CP 
violation induced by loop effects hardly leads to CP-violating Higgs couplings to fermions. Large 
CP-violating effects are generated by the CP-violating phases present already at tree level in 
the NMSSM Higgs sector. While the parameter points for non-vanishing values of are to 
a large extent excluded by the EDM bounds, those for non-trivial <p 2 values respect the EDM 
constraints. In this case, CP violation in the Higgs couplings to r leptons may be accessible 
exp er iment ally. 
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